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Abstract. We analyze a one dimensional quantum model with off-diagonal disorder, 
consisting of a sequence of potential energy barriers whose width is a random variable either 
uniformly or normally distributed. We investigate how the disorder and the energy distribution 
(due to a thermostat at room temperature) affect the resulting value of the transmission 
coefficient, and discuss the structure of the fluctuations of such coefficient at different length 
scales and the onset of different size-dependent regimes. Our analysis suggests an efficient 
way to detect tunneling resonances in a regime of off-diagonal disorder 
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1. Introduction 

Nonequilibrium thermodynamics [jTl is based on the vast separation between the space and 
time scales of the microscopic, mesoscopic and macroscopic physical realms. This separation, 
which characterizes the condition of local thermodynamic equilibrium |[D|2l[3l, cannot be 
realized in small systems, such as those concerning the modem bio- and nano-technologies. 
Hence, the standard statistical mechanical approach, based on the thermodynamic and 
hydrodynamic limits, should be extended to treat systems made up of relatively small numbers 
of microscopic components. In particular, one should account for the fluctuations of the 
physical properties of such systems, as they can be of the size of the average signals. In 
this work we face these issues considering a variant of the Anderson model H, which is the 
prototype of a disordered solid |l5l. In particular, we investigate the role of the microscopic 
disorder on the transmission coefficient of systems with a number N of potential barriers 
and in the presence of a thermal reservoir at a given temperature T. For this purpose, we 
consider one dimensional systems consisting of a sequence of barriers and wells of randomly 
chosen widths, under the constraint that the sum of the barrier widths and the total length 
of the system are fixed and do not change when grows. Consequently, our model enjoys 
a purely off-diagonal disorder [[6l |71 which affects the tunneling couplings among the wells, 
but not the energies of the bound states within the wells. This is not the case of the original 
tight-binding model, whose random fluctuations only concern the energy of a bound state 
[l8l|9l, introduced by Anderson in his pioneering paper [4] on localization effects in disordered 
solids. Furthermore, while in Anderson's model increasing the number of barriers corresponds 
to taking the large system limit, in our case it corresponds to distribute more finely the 
same amount of insulating material within the fixed length of the system. Therefore, our 
results differ substantially from those regarding the standard Anderson's model, and can be 
summarized as follows. 

• There are no localization effects for the equilibrium distribution of energies at room 
temperature: positive currents persist even in the large A^ limit; Furstenberg type 
theorems do not apply, flU . Mathematically, the reason is that the product of the first A^ 
random matrices yielding the transmission coefficient for a given choice of A^ barriers, 
changes, in order to preserve the length of the system and the sum of the barrier widths, 
when the A^ + 1-th matrix of the A^ + 1-th barrier is introduced. Physically, our models 
realize a kind of off-diagonal disorder which leads to delocalization. 

• The value of the transmission coefficient, averaged over the ensemble of disordered 
configurations, is close, for large A^, to the value corresponding to the ordered sequence 
of equally spaced barriers and wells, which is bounded away from zero. 

• The transmission coefficient increases non-linearly with the temperature parametrizing 
the distribution of the energies of the waves entering the system from the charge reservoir, 
which is also a heat bath. 

• There is a scale for A^, above which the (always positive) transmission coefficient no 
longer depends on the specific realization of the disorder, but still depends on A^, and 
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there is another scale above which it no longer depends on A^. We call "mesoscopic" the 
first, since it corresponds to nanoscopic samples, and we call "macroscopic" the latter 
scale, since it corresponds to macroscopic nanostructured materials. This means that all 
realizations of the disorder are equivalent in the N oo limit. 

• At room temperature, the probability distribution function (PDF) of the fluctuations of 
the transmission coefficient satisfies a principle of large deviations. Furthermore, the 
peak of this PDF corresponds to the transmission coefficient of the regular realizations. 

• Tunneling resonances for off-diagonal disordered systems in contact with an external 
bath can be predicted by investigating the corresponding ordered systems, which are 
amenable to an analytical treatment. 

2. The model 

We consider a model of a macroscopic semiconductor device (the bulk) consisting of an array 
of potential barriers and A^— 1 conducting regions (wells), in equilibrium with one electrode 
acting as an external thermostat at temperature T = 300/i (cf. Fig. [I]). The barriers have 
a constant height V{x) = V while their width is either uniformly or gaussian randomly 
distributed. For any N, the widths of the conducting regions take a constant value 5n- We 
denote by L the fixed total length of the sample, by L** the sum of the widths of all the barriers 
(i.e. the total length of the insulating region), and we also assume that the ratio (3 between 
insulating and conducting lengths takes the same value for all different sample realizations, 
so that 



where m is the mass of an electron. The boundary conditions prescribe Aq > for the 
amplitude of the plane wave entering from the left boundary and A^^^i = (no wave enters 
or is reflected from the right boundary). The barriers are delimited by a set of 2N points, 
denoted by Xq = 0, X2n-i = L in Fig. [H hereafter called nodes of discontinuity of 
the potential. The left boundary consists of a classical thermostat at temperature T, from 
which particles emerge at node in xq as plane waves, with energies distributed according to 
the Maxwell-Boltzmann distribution. Differently, no particles come from the vacuum on the 
right. 

Thus, denoting by Ui the i-th region, for £ G {0, 2, 2N}, the solutions of eq. take 
the form: 



L = {l + f3){N-l)S^ 
holds. This 1-D microscopic model obeys the steady state Schrodinger Equation: 



(2) 



(1) 





v42£e-^^ + Aa^+ie^^ ifxeUi and i is odd (i.e. for V{x) = V) 



—ikx 



if X eUi and i is even (i.e. for V{x) = 0) 



(3) 



with k = V2mE / h and z 



y/2m{V - E)/h. 
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Figure 1. 1-D Multiple-wells system, consisting of a sequence of 2N + 1 regions; N potential 
barriers, whose width is uniformly randomly distributed, {N — 1) conducting regions of 
constant width 6n, and the two boundary regions, on the left and on the right. The left 
boundary is modeled as a classical thermostat at temperature T. Thus, the mean energy of 
the plane waves entering from the left boundary is given by ^ksT. 



Following lfT3l . in each of the conducting regions, one may define the steady state 
currents as: 



2mi 



= jr{A,,)-jf{A,e^,) , (4) 

where the * denotes complex conjugation, jj^{A2e) = hk / m\A2i\^ denotes the current 
transmitted from the — l)-th barrier on the left (or, for £ = 0, from the thermostat located 
at the left boundary) and (^2^+1) = fi'k / m\A2i+i\'^ denotes the current reflected from the 
(£+l)-th barrier, cf. Fig. [T] Then, the application of the BenDaniel-Duke boundary conditions 
ifTTI . which require the continuity of ipt{x) and -^ipi{x) at the nodes, results in the constancy 
of the value je{x) across the wells and entails je{x) = i2N{x), for every even I. Equation dH), 
together with eq. leads to the following definition of the transmission coefficient S across 
the system: 

^ ^ " IfiAoT - IW ^ ^ 

which depends on the parameters of the model, such as the number of barriers N, the height of 
the barriers (which contributes to the parameter z) and the particular realization of sequence 
of barriers. In order to numerically compute the coefficient 5* as a function of the various 
parameters of the model, it proves convenient to rescale eq. ^ with respect to characteristic 
quantities, in order to rewrite it in a dimensionless form. For this purpose, let us introduce 
X = Lx, with L given by ©, ^ = EEt, iIj{x) = ^p{x){^/Ly^, V = VEt, with Et = ksT 
(i.e. twice the mean kinetic energy of the plane waves entering the bulk from the left side). 
Moreover, by introducing the scalar parameter a = h"^/ (2mL^i?T), one obtains the following 



expression for the dimensionless wavevectors: k = i/aV E and z = ^JayV — E. In the 
sequel we will refer to the dimensionless quantities and, to this aim, we may omit the tilde 
symbols, for sake of simplicity. The dimensionless version of eq. ©, then, attains the form: 

dx 



^ij{x) = a~\V - E)^Ij{x), xe[OA] (6) 
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which is the Schrodinger Equation (SE) which will be solved numerically with the 
aforementioned conditions at the nodes. 



3. The transfer matrix technique 

Let us describe our method of solution of the SE which follows Refs. [fTTl [T2ll and is 
referred to as the Transfer Matrix (TM) technique. Using eqs. the boundary conditions 
on the generic £-th node, with £ E {0, 1, 2N — 1} read as: 

where xe{x) = Y!IX K + 5Ar£/2 if i is even, and xe{x) = ^^ + ^n{£ - l)/2 if i 

is odd, and where Aj denotes the (dimensionless) random width of the i-th barrier. In matrix 
form, eqs. ^ can be written as: 



Mo(xo) 

M2(Xi) 



Ai 

A2 
^3 



A 




M4(a:2) ■ I = M5(X2) ■ I 7 1 (8) 



5 



M4Ar_2(x2Ar-l) ■ | ^ | = M4Ar_i (x2Ar-l' 




where, as long d& E <V , the 2x2 matrices of coefficients yi2t{xi) and M2£+i(x^) read 

M^eixe) = and M2,^,(x,) = ^^^^^ J (9) 



or 



yi2,{x,) = I ^^^^^ j and M2,+i(a;,) = -ike-^>^^^ ' ^^^^ 

if i is, respectively, even or odd. Assuming that the amplitude Aq of the incoming wave is 
known and that A^n+i = 0, because there is no reflection from the right boundary in the 2N- 
th region, then the linear system ([8]) corresponds to a set of AN equations in the AN variables 
{Ai, A^n}. Skipping the spatial dependence of the matrices to keep our notation simple, 
eq. ^ implies 
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^° J = Mq ^ ■ Ml ■ ^ ■ M3 ■ ... ■ M^^., ■ M4^_i ■ I 

= M.(^^J-j , (11) 

where we have also defined the 2x2 matrix M. Equation (fTTI) leads to 

Ao = MuA^^ (12) 

with Mil denoting the element on the top left corner of the matrix M. Relation (fT2l) allows 
us to rewrite eq. ([5]) in the form: 

A*,Ao -|Miip- ^''^ 
The results obtained with the TM method to the model described by eq. © are illustrated in 
SecH 



4. Numerical results 

We characterize the presence of disorder in the sequence of barriers letting p{X)dX be the 
probability distribution of the widths of a generic barrier, with i = {1, A^}, to take 
values in a range d\ centered on A. In particular, we used pseudo-random generators to 
investigate two relevant choices for p(A). The first is the uniform density, with support on 

the unit interval, while the other is the density p = \/^e~~, supported on R+, which is 
readily obtained from the gaussian density p = J\f{0, 1) by retaining only the positive values 
of the Aj's. Each of the two distributions is characterized by the corresponding mean (A) and 
variance a^. It proves useful to introduce, for both the distributions above, the realization 
mean and variance, defined, respectively, as 

TV 



Afi 

i=l 



In the large N limit, the random variable A b converges in probability to the mean (A) (Weak 
Law of Large Numbers), while the random variable converges with probability 1 to cr^. 
Since we use dimensionless variables in eq. we introduce the rescaled barrier width as: 

Xi = cXi, with c = — — . (14) 

{1 + (3)NXb 

I For a uniform density (A) = 0.5, ct^ = 1/ 12, whereas for the "gaussian" density defined above, with support 

on R+, (A) = 2/n, ct^ = 1 - 2/tt. 
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Therefore, for any given N and /3, the rescaled mean 



A 



B 



[1 + (3)N 



is no longer a random variable, and attains the same constant value independently of the 
chosen probability density. On the other hand, the rescaled realization variance 



^2 



N 



i=l 

remains a random variable which, for large N converges a. e. to = c^a^. We introduce the 
vector- valued random variable A at, defined by 

= {Ai,...,Ajv} , (15) 

which corresponds to a given realization of the sequence of barriers and will be referred to as 
a microscopic configuration. 

For given /3 and A^, one may, then, consider the collection = {A^-*, A^^'^-'} of A',, 
random realizations of the sequence of barriers which have been constructed numerically. 
Then, the average of a random observable O over the sample ^2 writes: 

(^)-=i^E^K'0 ^^^^ 

Among the possible configurations, the regular one 

A^ = {A^...,A^} , (17) 

which approximates the infinite superlattice of the literature on Bloch waves [fTOlfTTll . will be 
crucial also in our work. In our numerical simulations we set /3 = 0.1 and we investigated the 
behavior of the coefficient ^(A^, A^v, V, E; T) at a given temperature T, as a function of the 
number of barriers A^, of the microscopic configuration Ajy, of the (dimensionless) energy of 
the barriers V and of the (dimensionless) energy E. 

Figure [21 illustrates the behavior of S{N, An, V, E; T), S{N, A^, V, E; T) and of the 
average {S)n over an ensemble characterized by p(A) = 1 at T = 300K and for two different 
values of the wave energy E, one of which leads to a condition of resonant tunneling [fT2l|. 
The plots reveal that the values of S and of its average {S)q tend, for growing A^, to the same 
value, as intuitively expected. However, it is interesting to note that this value coincides with 
the most probable value of S{N, An,V, E;T), which is the value corresponding to the regular 
lattice configurations (flTl) . hereafter denoted as Sb{N, V, E; T). This holds for all values of 
V and of E. Let us now consider a thermostat located at the left boundary, so that the plane 
waves entering the bulk have an energy obeying a classical equilibrium distribution at a given 
temperature T. In the following plots we consider a one dimensional Maxwellian probability 
density 



f,,{E) = ^l/{TxE)er'^ (18) 
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Figure 2. Top panel: Behavior of S{N, An,V, E;T) vs. TV, for V = 20, E = 2.0 at 
T = 300K and for different microscopic configurations: the red curve corresponds to a single 
uniformly distributed configuration A^v, the purple curve represents the average (5*)^ over an 
ensemble of uniformly distributed configurations, while the dashed green curve corresponds 
to the realization ( fTTI l. Bottom panel: Same curves as those shown in the top panel, but here 
evaluated at E = 13.8, where a condition of tunneling resonance is realized. 

and we average over all energies to obtain 

POO 

SiN, A;v, V;T)= / 5(iV, A^v, V, E; T)f,,{E)dE (19) 

where the coefficient S, integrated in the r.h.s. of eq. (fT9l) . is computed according to eq. (fT3l) . 
The typical behavior of S{N, A^v, V] T) as a function of the number of barriers A^, for V = 
is illustrated in Fig. |3l which refers to three different microscopic configurations: the regular 
lattice A^ introduced in eq. (flTl) (red curve), a random uniform configuration (green curve), 
and a "gaussian" distribution (blue curve). 

The result illustrated in the top panel of Fig. [3] shows that, for large A^, the values 
of S{N , Kn,V]T = 300) corresponding to single random realizations fluctuate around 
Sb{N, V;T = 300) for both distributions. Similarly, the lower panel of Fig. |3] shows the 
behavior of ^(A^, A^v, V; T) vs. N for three different random uniform configurations. 
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Figure 3. Top panel: Behavior of S{N,An,V;T) for V ^ 20 and at T = SOOi^ 
corresponding to three different microscopic configurations: the regular lattice configuration 
(thick green dashed curve), a uniformly distributed configuration (red curve) and a 
normally distributed one (blue curve). Bottom panel.Behavioi of S{N, A]\[,V;T) vs. N, 
for V = 20, T ~ 30QK and for three different random uniformly distributed configurations 
Aj^\with^== 1,2,3. All curves tend to overlap with that corresponding to the regular lattice. 



We also investigated the dependence of S{N, An,V]T) on T. The top panel of Fig. |4] 
corroborates, in the limit of large N, the numerical results illustrated in the top panel of 
Fig. [3l Namely, the trend of the random values ^(A^, Aat, V; T) to approach Sb persists even 
at lower and higher temperatures. This is one of the prominent features of our model. The 
lower panel of Fig. |4] shows that the variation of Sb, over a wide range of temperatures, can 
be suitably fitted by a quadratic curve, 0(T) = aT^ + 6T + c. 

Figure [3] anticipates two further crucial aspects which will be addressed in more detail 
below. The first concerns the magnitude of the fluctuations of the values of S{N, Ajy, V; T), 
(iV, V) = ^{{S- {S)ny)n, which decreases with N. This decay of the size of the 



a 



fluctuations allows us to identify a "mesoscopic" scale Nmeso- The number Nmeso depends on 
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Figure 4. Top panel: Values for the observable S'(A'^, Aat, T) at different temperatures, 
with V = 20. Shown are the numerical results for uniformly distributed configurations 
(straight lines) and for the corresponding regular lattice configurations (thick dashed 
curves). Bottom panel: Fitting of the numerical data of Sb{N, An, V), calculated at = 300 
and with V ~ 20, at different temperatures with the curve 4){T) = aT^ + bT + c, with 
a = -2.022 ± 1.891 • IQ-'^, b = 0.002 ± 9.672 • lO^^ and c = -0.237 ± 0.011. 



V and the coefficient S{N, Ajv, V; T) depends only weakly on the microscopic configuration, 
if > Nmeso- The second scale Nmacro > Nmeso is such that S{N, A]\f,V;T) depends 
neither on the configuration nor on the number of barriers if > Nmacro- In other words, S 
attains a value which is, practically, that concerning infinitely finely structured samples at a 
finite fixed length. Note that Nmacro depends also on V. 

Next, in Fig. [51 we plotted the behavior of {S)n{N, V] T) with A^, for V = 20, and 
we compare it with the value of Sb{N, V; T) for regular lattices. The result is consistent 
with those of Figs. [2]and|3l in that it shows the regime corresponding to A^ ^ 1, where the 
curve of {S)n{N, V; T) varies very slowly with A^ and approaches the value Sb{N, V; T). 
This convergence process is a collective effect, in that it does not depend on the underlying 
disorder and is characterized by the delocalization of the wave function typical described by 
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Figure 5. Behavior of {S)n vs. N, for V = 2Q,T = 300A' and for different values of Nr. 

the Bloch waves theory [jlOil for regular lattices. The plot in Fig. [5] reveals that, given an 
ensemble Vt of random uniform configurations, one has: 

{S)n{N,V;T)c^SB{N,V-T) for large iV . (20) 

Indeed Fig. [5] shows that N > 200 suffices to reach a good accuracy even with small samples 

n. 

Moreover, we verified that the energies of the incoming wave, at a given temperature, 
lie in the conducting band of the infinite periodic chain of barriers, which implies that 
SBiN,V;T) > 0. A few comments are in order, here. In the first place, the absence of 
localization can be traced back to the fixed finite amount of insulating material, which we 
have even in the N ^ oo limit, because L is fixed. As a consequence, incoming waves may, 
at most, be damped by a finite factor, except, perhaps, for a negligible set of energies which 
we have not observed. This distinguishes our model from the tight-binding model, which is 
more extensively investigated in the specialized literature, and also prevents the application of 
the Furstenberg's theorem Q. Indeed, introducing the (N + l)-th barrier in one our systems 
amounts to a complete rearrangement of the previous N barriers. Mathematically, this means 
that the product of the first N random matrices is replaced by a new product. Differently, the 
case of ergodic-like theorems, such as Furstenberg's one, applies to products of N random 
matrices which are multiplied by the (A^ + l)-th matrix. Moreover, the decrease of the size 
of fluctuations with A^, which will be explored in more detailed below, can be regarded as a 
phenomenon of self-averaging of the observable S, which is quite common [[5]|. However, our 
results, further supported by the analysis of the PDF of the transmission coefficient. Fig. 16 17 1 
below, show that the average over the disorder of the random values ^(A^, A^v, V; T) attains 
the specific value Sb{N, V; T) pertaining to the regular configuration. In other words, the 
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transmission coefficient self-averages around a value {S)q{N,V;T) which is close to Sb 
even for small and tends to Sb in the N ^ oo limit. As shown in Fig. Ul given a sample Q 
of uniform realizations, Sb{N, V; T) corresponds to the most probable value of the random 
variable S in the sample, which, when N grows, tends also to the mean {S)q{N,V;T). 
Furthermore, being obtained from the dimensionless SE this result is independent of 
the length of the system (as long as L is fixed) and does not result in localization effects when 
L is large, because Sb does not vanish. 

Let us now investigate, more accurately, the structure of the fluctuations, in the sample Q 
of random uniform configurations at temperature T = 300Ji. Denote by pj^(S)dS the 
probability that S falls in the interval dS, so that: 

{S)n = [ Spj,iS)dS . (21) 



The numerical results presented so far on the relation between {S)n{N,V;T) and 
Sb{N, V; T) and their mutual relation, eq. (|20|) . as well as on the decrease of the fluctuations 
size with growing N, indicate that pn{S) peaks more and more around the reference value 
Sb{N, V] T). To show this, we numerically calculated the quantity pn{S) from eq. (fT6l) and 
we plotted in Fig. [6] the resulting curves for different values of and for V = 20. The 
maxima of the PDF in Figl6]are located at S = ^^(A^, V; T), cf. also Fig. Ul and tend, for 
large A^, to the mean values {S)q{N, V; T). One further realizes that pAr(S') obeys a large 
deviation principle with a given rate functional (j^, in the sense that the limit 

N^oo iV Af— i-oo 

exists. Figure |7] shows, for the range in which we have good statistics, that the large deviation 
rate functional is apparently smooth and strictly convex. It is worth pointing out, again, that 
the A^ — 7- oo limit is not achieved in the standard fashion of products of random matrices. 

It is also interesting to note, in Fig. Ul that the locus of the minima, i.e. of highest 
probability density, of all the curves (NiS) is represented by SsiN, V; T). Moreover, Fig. [8] 
reveals that, for small deviations from ^^(A^, V; T), the rate functional CNiS) is quadratic, 
as expected where the central limit theorem applies. Following the approach outlined in Ref. 
ifTSl . if a large deviation functional ({S) exists, for any finite A^ one may write: 

C{S) = ^\ogp^iS) + zr,iS) (23) 

where is a correction term expected to be of order 0(;^) as long as correlations decay 
sufficiently fast, cf. lfT6l ITTl . This correction term may be further split as zn{S) = 
(p{S) /N + zpf{S). Thus, if we consider two different values of A^, say Ni and N2, one may 
write ifTSll : 

^^^^ = N^-N2 /'^'^^ ^ ^ ^'^"^^ " ^'^""'^ ^^"^^ 

where the remaining correction terms zj^- are expected to give a smaller contribution to than 
zn in eq. (|23l) does. We numerically verified the validity of the finite A^ estimate of the large 
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Figure 6. Probability densities pn{S) for different values of N and for V = 20. 




Figure 7. Rate functional associated to pn{S) for different values of N and for V ~ 20. 
The curves Ctv move downwards for growing N, so that, as expected, in the iV oo limit, 
C{S) intersects the horizontal axis only in 5* = Sb- 
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Figure 8. Rate functional Cn{S) associated to pn{S). Blue points: results of the numerical 
simulation with N = 100 and V = 20. Green dashed line: Fitting of numerical data 
with the parabola 0(5* - Sb{N)Y + b, with parameters a = 2635.5 ± 82.96 and h = 
-5.655- 10-2 ±3. 38 - 10-4. 
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Figure 9. Plot of the function ivT^ log (^) with A^i = 125 and N2 
{10,25,50,75,100}. 
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Figure 10. Top panel: Fit of Opj Sm (red curve), evaluated for V = 20, with the curve 
f[N) = a X N'K with a = 3.798383 x 10"^ ± 4 797 x 10"'*. The plot confirms the 
asymptotic behavior of (Jp/Sm, given by eq. ( |25] l. Bottom panel: Behavior of <Tp{N, V) vs. 
A'^, for V = {25, 30, 35}. The thick orange dashed line denotes the value cr* below which we 
may conventially assume that the size of fluctuations is practically negligible. 



N] , for iV > 1. (25) 



deviation functional C('S'), given by eq. (|24|) . observing the rapid collapse of all curves of Fig. 
|7]into a single curve, in Fig. [9l 

From the validity of a principle of large deviations and of the central limit theorem, one 
expects the following asymptotic behavior for the fluctuations of the transmission coefficient: 

We numerically checked that this is the case by evaluating the ratio in eq. ( |25] ). cf. the top 
panel of Fig. \T0[ for V = 20. The bottom panel of Fig. \T0\ further illustrates the decay 
of Nmeso with V: taking N„ such that ap{N,V) < <y*p = 1-4 x lO'^ VA^ > Nrr,eso, 
where a* is considered small, Nmeso{V) is found to rapidly decrease with V: we obtain 
Nmeso = 100, 18, 8 for, respectively, V = 25, 30, 35. Then, Nmeso tends to when the barrier 
height diverges: in this limit the transmission coefficient vanishes for any configuration A^v, 
hence the fluctuations are suppressed. 
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5. Concluding remarks 

In this work we investigated a quantum mechanical one-dimensional system, in equilibrium 
with a classical heat bath. The role played by the disorder on the transmission coefficient S 
has been investigated. Our numerical investigation sheds light on the existence of appropriate 
scales, respectively delimited by Nmeso and N macro, which are not as widely separated as 
in thermodynamic systems. The main result of our work suggests that, in presence of an 
external heat bath and off-diagonal disorder, the wave function is delocalized over all the 
system length. This absence of localization reflects the long range correlations which extend 
over the whole length of our systems and prevent the low tunneling limit. Therefore, the 
effective tri-diagonal Hamiltonian for given bound states in a potential well, which is common 
to studies of the Anderson model, does not pertain to our systems and should in case be 
replaced by a full matrix. The presence of correlations is known to lead to extended states 
in models enjoying off-diagonal disorder, cf Ref. JH dH. Furthermore, our procedure for 
the N ^ oo limit makes the Furstenberg theorem inapplicable. The variable S is self- 
averaging for growing A^, and tends to the most probable value. Interestingly, this value is Sb, 
which corresponds to a microscopic ordered array of barriers and wells. This observation has 
technological implications, as it suggests the possibility of predicting tunneling resonances 
for off-diagonal disordered systems in contact with an external bath, by investigating the 
corresponding ordered system, which is amenable to an analytical treatment. Moreover, we 
have numerically observed the validity of a principle of large deviations for the probability 
density pn{S), extending the results presented in Ref. [7], where, finally, we verified that the 
most probable value of the transmission coefficient corresponds to the regular configurations. 
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